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ABSTRACT 

We present Spitzer IRS spectra of four carbon stars located in the Galactic Halo and 
the thick disc. The spectra display typical features of carbon stars with SiC dust emission 
and C2H2 molecular absorption. Dust radiative transfer models and infrared colors enable us 
to determine the dust production rates for these stars whilst prior CO measurements yield 
expansion velocities and total mass-loss rates. 

The gas properties (low expansion velocities (around 7 km/s) and strong C2H2 molecular 
absorption bands) are consistent with the stars being metal-poor. However the dust content 
of these stars (strong SiC emission bands) is very similar to what is observed in metal-rich 
carbon stars. The strong SiC emission may indicate that the carbon stars derive from a metal- 
rich population, or that these AGB stars produce silicon. 

The origin of the halo carbon stars is not known. They may be extrinsinc halo stars be- 
longing to the halo population, they may have been accreted from a satellite galaxy such as the 
Sagittarius Dwarf Spheroidal Galaxy, or they may be escapees from the galactic disk. If the 
stars are intrinsically metal-rich, an origin in the disc would be most likely. If an a-element 
enhancement can be confirmed, it would argue for an origin in the halo (which is known to be 
a-enhanced) or a Galactic satellite. 

Key words: circumstellar matter - infrared: stars — carbon stars — AGB stars — stars: mass 
loss 



1 INTRODUCTION 

Low- and intermediate-mass stars (LEVIS, 0.8-8 M© on the main 
sequence) end their lives with an intense mass-loss episode on 
the asymptotic giant branch (AGB). The mass loss from AGB 
stars is one of the main sources of enrichment of the interstellar 
medium (ISM) with newly synthesized elements. In our Galaxy, 
mass loss f rom AGB st ars produces well over half of the mea- 
sured dust jGehrzlll989l) . The AGB also dominates the measur- 
able dust production in the LMC a nd SMC (Matsuur a'et alj|2009l ; 
ISrinivasan et all2009l ; lBover et al.ll2012t) . 

The mass-loss mechanism is not fully understood, but a two- 
step mechanism is widely accepted as plausible (e.g., Winters et al. 
2000). Pulsation from the star extends the atmosphere, where the 
material is dense and cool enough to form dust. Radiation pressure 
drives the mass loss, accelerating the grains while the gas is car- 
ried along by friction. Dust formation is the key step in the mass- 
loss process. Its composition depends on the C/O ratio. If C/O >1, 
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i.e. the star is carbon-rich, then most of the oxygen is trapped in 
volatile CO molecules, which have high dissociation energy. The 
dust forms from the remaining refractory elements, and consists of 
amorphous carbon and SiC. For C/O < 1, the dust will be oxygen- 
rich, consisting of silicates and aluminium oxides.. 

Understanding the effect of metallicity on the mass loss is of 
prime importance to determining the yields of LIMS in different 
galaxies and understanding the formation of dust in the early Uni- 
verse. We thus carried out several surveys of mass-losing AGB 
stars in the Local Gro up with the Spitzer Infrared Spectrograph 
(IRS iHouck et ail 2004» targeting the Large Magellanic Clou d 
(LMC lZiilstra et al.ll200dlMatsuura et al. 2006; S loan et al.ll2008h. 
Smal l Magellanic Cloud (SMC ISloan et a l. 2006; Lagade c et al.l 
2007; Sloan et al. 2008), and several dwarf spheroidal galax- 
ies jMatsuuraetalJl2007l : lLagadec etal]|2009l ; I Sloan et al.l [20091 

The results show that dust-production rates from oxygen- 
rich AGB stars are lower in more metal-poor environments, 
but for carbon stars, m e tallicity shows no strong influence 
fcroenewegen et all |2007| ; ISloan et al] l2008h . Subtantial dust- 
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production rates are observed for carbon stars with metalliciti es as 
low as ~0.1 Z Q dLagadec et alj2008l ; lsioan et al J2009ll2012h . 

Most of the material ejected via the mass-loss process is in 
the gaseous form. The gas mass in the envelope of an AGB star 
is typically two orders of magnitude higher than the dust mass to 
which our infrared surveys are sensitive. To fully characterize the 
mass loss from AGB stars and its dependence on metallicity, one 
needs to study the gas mass loss in environments with different 
metallicities. 

Carbon stars have recently been disco vered in the Galactic 
Halo. These are thought to be metal-poor (Totten & Irwin 1998; 
lMauronetai]|2004 120051. 120071) . iGroenewegen et al.l dl997|) ob- 
served the CO emission from a carbon star in the Halo. We recently 
observed a sample of carbon Halo stars in the CO J = 3 — > 2 
transition and found that their expansion velocities were low com- 
pared to AGB stars in the Galactic disc with similar infrared colors 
i lLagadec et all2O10h . 

Here, we report new Spitzer IRS observations of four of these 
carbon stars. We place these observations in the context of our CO 
observations and additional IRS observations of three other likely 
carbon Halo stars. 



2 THE SAMPLE 



Lagad ec"et al.l ([2010) observed six carbon stars in the Halo us- 
ing the James Clerk Maxwell Telescope (JCMT) to study the CO 
J= 3 — > 2 transition at 345 GHz. Those s tars were select e d from 
catalogs of carbon stars i n the Halo by iTotten & Irwinl J 1998h . 
iMauron et all fe(X)4 12001 l2007n . and lMauroij d2008h . Each target 
had to be a bright source as observed by the Infrared Astronomical 
Satellite (IRAS) and have &J — K color of ~3^1. The near-infrared 
color requirement ensured that all of the dust shells had similar op- 
tical depths. 

We selected five of the six JCMT targets as part of an IRS pro- 
gram for the final months of the cryogenic portion of the Spitzer 
mission. We obtained spectra for four of our sample before the 
cryogens were exhausted, with the fourth spectrum coming on 2009 
May 15, the day that Spitzer ran out of liquid helium. In fact, this 
spectrum (of IRAS 08427) was the penultimate spectrum ever ob- 
tained by the IRS. 

In this paper we also examine three additional carbon stars 
likely to be in the Halo. lSloan et al. ||201(3) published the IRS spec- 
trum of Lynga 7 VI as part of their sample of AGB stars in globular 
clusters. Carbon stars are rare in globular clusters, making it possi- 
ble that this star is in the foreground or background of the cluster. 
However, all estimates of its distance coincided with the position 
of the cluster. Lynga 7 VI could be either a carbon Halo star in the 
immediate vic i nity o f Lynga 7, or it could actually be a member. 
lLagadec et all (2009) published the spectra of several carbon stars 
in the Sagittarius Dwarf Spheroidal galaxy (Sgr dSph) (previously 
known as the Sagittarius Dwarf Elliptical Galaxy (SagDEG)), in- 
cluding IRAS 18384-3310 and IRAS 19074-3233 (herafter IRAS 
18383 and IRAS 19074). The radial velocities of these two sources, 
as determined from optical spectra, showed that they were not 
members of Sgr dSph, and distance estimates placed them at least 
10 kpc in the foreground. 

Table[T]gives the observational properties of the nine probable 
carbon Halo stars observed with either JCMT or the IRS. For the 
six Halo stars with JCMT observations, the classifications of their 
membership in the th ick disk, Halo, or Sagittarius Stream are from 
lLagadec et all (2010), who used their radial velocities, distances, 



and galactic coordinates as indicators. The last two targets in the 
table, IRAS 18384 and IRAS 19074, are likely to be in the Halo, 
based on their radial velocities (—23 and —73 km/s, respectively). 
Lynga 7 VI is only 400 pc from the Galactic plane, consistent with 
membership of the old disc or thick disc. We will revisit these as- 
signments in Sec. 5.1 below. 

All IRS observations used both the Short-Low (SL) and Long- 
Low (LL) modules to cover the wavelength range 5-37 fim at res- 
olutions of ~100. Before extracting spectra from the IRS images, 
we (1) differenced them to remove background emission and ad- 
ditive spikes and divots, and (2) replaced pixels suspected of mul- 
tiplicative spikes and divots with values calculated from neighbor- 
ing pixels. We used the standard tapered-column extraction algo- 
rithm available with CUPID to obtain the spectra from the im- 
ages Qwhen combining spectra from the two nod positions in each 
aperture, we flagged and ignored bad data which appeared in only 
one of the nods. The spectra were then corrected with scalar mul- 
tiplication to remove discontinuities between spectral segments. 
These corrections were generally between 5 and 10%, but for IRAS 
04188, the two SL segments were raised by 21 and 13% (for the 5- 
7.5 and 7.5-14 fj,m segments, respectively), probably due to a slight 
mispointing. The spectral flux standards were HR 6348 (K0 III) in 
SL and HR 6348 and HD 173511 (K5 III) in LL. For more infor- 
mation, we refer the reader to the more detailed descriptio n of the 
observing method and data reduction by Slo an et alj d2012» Fl 



3 DESCRIPTION OF THE SPECTRA 

Fig.Q] displays the four infrared spectra obtained in our observ- 
ing p rogram, along with photometry from 2MASS dCutri et all 
2003) and the W ide-field Infrared Survey Experiment (WISE 
Wright et al] d2010h ). The spectra confirm all four sources as dusty 
carbon stars. Each shows a prominent SiC dust emission feature at 
11.3 fim, as well as absorption bands from acetylene (C2H2) at 7.5 
and 13.7 /im. The absorption band only partially covered by the 
IRS at 5 fim is due to CO and is often strong in carbon stars. The 
continuum between these various features and bands is primarily 
due to emission from amorphous carbon grains. Some sources also 
display a broad 30 fim feature. 

To analyze these spect ra, we apply the M anch ester Method, 
which was introduced by Sloan et alj d2006h and IZiilstra et al.l 
(2006). This method determines two colors from the spectra in four 
narrow bands. The [6.4] — [9.3] color is a good estimate of the opti- 
cal depth of the warm dust, while the [16.5] — [21.5] color provides 
an estimate of the overall dust temperature. For all the observed 
stars, we measured the strength of the spectral features from SiC 
dust at 1 1 .3 and acetylene gas at 7.5 and 13.7 /im using line seg- 
ments to estimate the continuum. In each case, we defined the con- 
tinuum using wavelength intervals to either side of the feature. We 
also measured the strength of the MgS dust emission at ~30 /im, 
but because the red edge of the feature is outside the IRS spectral 



1 CUPID is the Customizable User Pipeline for IRS Data, and it is available 
from the website of the Spitzer Science Center. 

2 The only difference between their methodology and ours is that the 
brightness of our sample did not require the use of the optimal extraction 
algorithm introduced bv lLebouteiller et alj ( 2010). 

3 We will refer to this feature as the "MgS" feature, although the reader 
should be aware of the continuing debate about whether or not MgS pro- 
duces the 30-/xm feature. See, for example, Hony et al. (2002), Zhang et al. 
(2009) and Lombaert et al. (2012). 
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Table 1. The obse rved sample of carbon stars in the Galactic Halo. The near-infrared photometry is from 2MASS, except for Lynga 7 VI, which is from 
ISloanetaljfioich . 



Name 




RA 


Dec 


Gal. coordinates 


K s 


J-K B 


Fl2 


Distance 


Scale height 


Possible 


Observations 






(J2000) 


1" 




(mag) 


(mag) 


(Jy) 


(kpc) 


(kpc) 


population 




IRAS 04188- 


-0122 


04 21 27.25 


+01 29 13.4 


192.2 


-32.0 


6.42 


3.28 


3.37 


6.5 


-3.4 


Thick Disc 


JCMT/IRS 


IRAS 08427+0338 


08 45 22.27 


+03 27 11.2 


223.5 


+26.8 


6.26 


3.41 


6.50 


5.5 


2.5 


Thick Disc 


JCMT/IRS 


IRAS 11308- 


-1020 


11 33 24.57 


-10 36 58.6 


273.7 


+47.8 


4.57 


3.83 


57.37 


2.5 


1.9 


Thick Disc 


JCMT 


IRAS 12560- 


-1656 


12 58 33.50 


+ 16 40 12.0 


312.3 


+79.4 


7.82 


3.48 


0.77 


12.0 


11.8 


Sgr Stream 


JCMT 


Lynga 7 VI 




16 11 02.05 


-55 19 13.5 


328.8 


-02.8 


7.25 


4.10 


2.99 


7.3 


-0.4 


Thick Disc 


IRS 


IRAS 16339- 


-0317 


16 36 31.5 


-03 23 35.0 


12.7 


+27.8 


6.10 


3.91 


14.57 


4.9 


2.3 


Halo 


JCMT/IRS 


IRAS 18120- 


-4530 


18 13 29.6 


+45 31 17.0 


73.1 


+25.3 


6.71 


3.81 


7.86 


6.7 


2.9 


Halo 


JCMT/IRS 


IRAS 18384- 


-3310 


18 41 43.50 


-33 07 16.6 


1.9 


-12.6 


8.58 


3.00 


0.80 


15.0 


-3.3 


Halo 


IRS 


IRAS 19074- 


-3233 


19 10 39.87 


-32 28 37.3 


5.0 


-17.9 


8.35 


3.61 


1.73 


14.5 


-4.5 


Halo 


IRS 



range, we used a blackbody with the temperature derived from the 
[16.5] — [21.5] color to extrapolate the continuum under this dust 
feature. The strength of the dust features are reported as the ratio 
between the integrated flux in the feature and the continuum un- 
derneath it. For the SiC feature, we also report its integrated flux, 
scaled to a distance of 10 pc. For the acetylene bands, we report the 
equivalent widths. Table|2]presents all of these results. 

To compare the strength of the observed features with ones in 
environments with different metallicities, we used a control sample 
of stars observed with the IRS and the SWS on the Infrared Space 
Observatory (ISO. This sample includes spectra from AGB stars in 
the Galaxy, LMC, SMC, Sgr dSph, and Fornax (see the references 
in Sec. 1). In all of the plots presented here, we use open symbols 
for Galactic stars. The symbols from largest to smallest represent 
stars in the Halo, the Thick Disc, and the Galactic SWS sample, 
respectively. Uncertainties were determined when coadding spec- 
tra extracted from individual images and propagated through all 
subsequent measurements. The uncertainties in the strengths of the 
feature in the IRS spectra include the effect of the uncertainty in 
the fitted continuum. 

Fig.[2]shows the strength of the SiC and MgS dust features as a 
function of the [6.4] — [9.3] color and the temperature derived from 
the [16.5] — [21.5] color. The strength of the MgS feature is low for 
all the observed stars, most likely due to the relatively high dust 
temperatures of the Halo stars. MgS start s to form around 600 K 
and completes its formation around 300 K dNuth et al.ll9 85). Fig.|2] 
shows that the dependence of the SiC feature on the dust tempera- 
ture is less clear. SiC tends to be weaker at high temperature (be- 
tween 1000 and 1500 K) than at low temperature (between 500 and 
1000 K). 

The upper left panel in Fig.[2] shows that the Galactic carbon 
stars follow one sequence in SiC/continuum strength as a func- 
tion of [6.4] -[9.3] col or, while the stars in the LMC, SMC, and 
Fornax follow another. lLagadec et al.l d2007t) explained these two 
sequences as a consequence of differences in the condensation of 
amorphous carbon dust at different metallicities. The Halo carbon 
stars in our sample follow the same sequence as the Galactic car- 
bon stars, indicating a similar dust formation sequence and perhaps 
a similar metallicity as well. 

The lower left panel in Fig.[2]shows the integrated fluxes of the 
SiC features, normalized to a common distance of 10 pc. While the 
separation between the two sequences is much smaller compared 
to the upper left panel, it is still apparent. Again, the Halo stars in 
our sample tend to follow the Galactic sequence. This plot removes 



the influence of the amorphous carbon continuum on the measured 
strength of the SiC feature and shows that in absolute terms, the 
feature is weaker in more metal-poor carbon stars. 

Fig.[3] displays the equivalent width of the 7.5 and 13.7 /jm 
C2H2 bands as a function of [6.4] — [9.3] color. In our small sam- 
ple, the Halo stars have stronger 7.5 acetylene bands than the 
stars assigned to the thick disc. To compare these stars to the other 
populations, it is necessary to limit the comparison to stars with 
similar [6.4] — [9.3] colors, since these have similar column densi- 
ties of amorphous carbon dust around them. The two Halo stars 
show band strengths much like those in metal-poor environments 
like the SMC or Fornax, while the two thick disc stars are more 
similar to carbon stars in the LMC and the Galaxy. The compar- 
isons of the 13.7 /im feature are less instructive, because we only 
measure the narrow Q branch of the feature, which is saturated for 
high column densities. 

Most of the spectroscopic studies based on IRS data ref- 
erenced in this paper have noted the general tendency for the 
7.5 ami acetylene ban d to strengthen in more metal-poor samples. 
iMatsuura et alj |2007) suggest that this trend is the result of higher 
C/O ratios in the outflows from the stars, but ISloanet "all ( 120121) 
have suggested it could result from less efficient dust formation. 
If high-temperature condensates such as TiC act as seeds for pro- 
duction of amorphous carbon grains from acetylene gas, the lower 
abundance of heavy elements like Ti in metal-poor carbon stars 
would lead to an excess of carbon-rich gas to carbon-rich dust. 



4 DUST-PRODUCTION RATES AND GAS MASS-LOSS 
RATES 

4.1 Dust-production rates 

To determine the dust-production rate (DPR, or D) from ou r car- 
bon stars, we follow the approach of lLagadec et af] (2010), who 
use two techniq ues. The first method uses near- and mid-infrare d 
flux of the stars jWhitelock et alJl994ll2006];|Lag"adec et alj2008h . 
Whitelock et al. I d2006l) give the following relation: 



log(M t , 



= -7.668 + 0.7305(_K" - [12]) 
-5.398 x IQ~ 2 {K- [12]) 2 
+1.343 x IQ- 3 {K- [12]) 3 , 



(1) 



where K is the 2MASS K s magnitude, and [12] is the IRAS 12 
magnitude, assuming a zero-magnitude flux of 28.3 Jy. We estimate 
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Figure 1. Spitzer IRS spectra of the four Halo carbon stars we observed. Filled squares and circles are 2MASS JHK S and WISE photometry. The dashed line 
represents the best fitting radiative transfer model. 



Table 2. Spectral measurements for our sample: The two colors, the equivalent widths (EW)of the acetylene bands, the SiC dust emission strength (normalized 
to the continuum and normalized to 10 pc), the MgS dust emission strength (normalized to the continuum), and the dust temperature, derived from the 
[16.5] — [21.5] color. 









EW (7.5 H 


EW(13.7 /im) 




SiC flux at 10 pc 




T 


target 


[6.4] -[9.3] 


[16.5] — [21.5] 


(/an) 


(fim) 


SiC/cont. 


(10" s Wm~ 2 ) 


MgS/cont. 


(K) 


IRAS 04188 


0.417±0.009 


0.212±0.011 


0.118±0.011 


0.054±0.008 


0.239±0.007 


1.30±0.04 


0.096±0.014 


639 ± 29 


IRAS 08427 


0.557±0.007 


0.176±0.012 


0.073±0.005 


0.032±0.003 


0.304±0.005 


2.79±0.05 


0.030±0.015 


754 ± 45 


Lyngaa 7 VI 


0.549±0.007 


0.197±0.013 


0.011±0.004 


0.030±0.003 


0.383±0.009 


4.91±0.11 


0.396±0.016 


683 ± 39 


IRAS 16339 


0.494±0.010 


0.163±0.015 


0.162±0.005 


0.061 ±0.005 


0.345±0.006 


4.21±0.08 


0.137±0.019 


813 ± 65 


IRAS 18120 


0.771 ±0.007 


0.170±0.014 


0.226±0.006 


0.030±0.003 


0.296±0.004 


3.50±0.04 


0.002±0.017 


782 ± 54 


IRAS 18384 


0.382±0.011 


0.273±0.010 


0.141±0.003 


0.035±0.003 


0.267±0.005 


0.69±0.01 


0.655±0.013 


510 ± 16 


IRAS 19074 


0.644±0.005 


0.097±0.016 


0.088±0.003 


0.030±0.002 


0.409±0.006 


1.22±0.02 


0.062±0.020 


1335 ±188 
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Figure 2. Left: Strength of the SiC dust emission feature as a function of the [6.4] — [9.3] color, with the top panel showing the strength normalized to the 
underlying continuum and the bottom panel showing the integrated flux normalized to a distance of 10 pc. Right: Strength of the dust features as a function of 
dust temperature (as estimated from the [16.5] — [21.5] color). The large open circles are the seven stars from the present study (the largest ones represent Halo 
stars, while the intermediate ones represent thick disc stars). The small open circles and other symbols are from the comparison samples, as shown in the inset 
keys. 



the DPR from the total mass-loss rate by dividing by the assumed 
gas-to -dust ratio of 200. The relation obtained by IWhitelock et all 
(2006) assumes an expansion velocity of 19 km/s. We scale our 
results using the expansion velocity as measured by the JCMT. 

The second method uses radiative transfer models. Here, we 
present new models using the Spit zer spectra as a n additional con- 
straint that was not available when lLagadec et alj ( 120101 ) generated 
the earlier models . We used the radi ative transfer code D USTY 
dlvezic etal.lll999h . as described by lLagadec et all ^20 10). This 
code solves the one-dimensional radiative transport problem in a 
dusty environment. All of the models assume that the irradiation 
comes from a point source (the central star) at the center of a spher- 
ical dusty envelope. The circumstellar envelope is filled with ma- 



terial from a radiatively driven wind. All of the stars are carbon- 
rich, and the dust consists of amorphous carbon and SiC. O ptical 
properties for th ese dust grains are taken from lHannerl fl983) and 
Pegourie(1988), respectively. The grain size distribution is taken as 
a typical MRN distribution, with a grain size a varying from 0.0005 
to 0.25 fim d istributed accordin g to a power law with n(a)oca -9 , 
where q=3.5 iMathi s et alll977h . We are aware that this grain size 
distribution was estimated for the interstellar medium, where dust 
grains are thought to be significantly processed, but this is as good 
an assumption as any given the lack of information about dust grain 
sizes in these stars. The outer radius of the dust shell was set to 10 3 
times the inner radius; this parameter has a negligible effect on our 
models. 
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Figure 3. Left: Equivalent width of the 7.5 fim C2H2 feature as a function of the [6.4] — [9.3] color. Right: Equivalent width of the 13.7 fim C2H2 feature as 
a function of the [6.4] — [9.3] color. The symbols are as defined in Fig.[3] 



To model the emission from the central st ar, we used a hy- 
drostatic model including molecular opacities dLoidl et all 1 200 ll ; 
iGroenewegen et al]|2007t) . Our aim was to fit the spectral energy 
distribution, defined by the 2MASS photometry (J, H and K s ) and 
our Spitzer data to estimate the DPRs. Our models included only 
dust, so we do not attempt to fit the molecular absorption bands 
seen in the Spitzer spectra. We fixed the dust temperature at the in- 
ner radius to be 1200 K and the effective temperature of the central 
star, T e fl = 2800 K, unless a satisfactory fit could not be obtained 
with these parameters. 

DUSTY gives total (gas+dust) mass-loss rates assuming a gas- 
to-dust ratio of 200. The DPR can be obtained by dividing by this 
value, assuming that the gas and dust have the same expansion ve- 
locities. The expansion velocity is an output for the DUSTY mod- 
els, and we scaled our results using the expansion velocities mea- 
sured with JCMT and assumed it to be 10 km.s -1 when no mea- 
surements were available. 

Table[3] gives the results from the fitted DUSTY models. Note 
that where available, t he DPRs determined f rom the K s — [12] color 
are simply taken from lLagadec et all J20ld> . 



4.2 Gas mass-loss rates 

Spectroscopic observations of the CO lines are our best means of 
determining the gas mass-loss rates from AG B stars. CO is , after 
H2, the most abundant circumstellar molecule d01ofssoiJl997h and 
can be used easily as a probe of the circumstellar medium. Since the 
1980s, numerous works, combining such observations and radiative 
transfer models, have attempted to find a formul a linking the CO 
line p roperties and the gas mass-loss rates (e.g. iKnapp & Morris! 
1985). From this point on, we will equate the total mass-loss rate 
with t hat determined from the gas. 

iRamstedt et alj ((2008) presented a reliable and up-to-date for- 
mula to determine the total mass-loss rate from CO observations: 

M = Sj(Ico0 2 b D 2 ) a >v b e >f-% , (2) 



where Ico, D, V e and fco are the CO line intensity (in Kkm s _1 ), 
the distance (in kpc), the expansion velocity (in kms" 1 ) and the 
CO abundance, respectively. The parameters Sj, aj, bj and cj 
depend on the CO transition observed. For the CO J = 3— > 2 line, 
they are 3.8x 10 -11 , 0.91, 0.39 and 0.45, respectively. 6 mb is the 
full width at half maximum (FWHM) of the main beam of the tele- 
scope, which is 14 arcsec at JCMT. We adopted the standard so- 
lar abundance of 8.78 x 10~ 4 for fco (Asplund et al. (2009)). This 
value could be lower as, in metal-poor environments, less oxygen 
is present to form CO, making the CO abundance with respect to 
H2 lower. That would make the total mass-loss rates higher than 
what we have determined below. Thus, the total mass-loss rates we 
estimate are lower limits. 

Table [4] shows the properties of the CO lines observed by 
lLagadec et al] (2010) and the total mass-loss rates they derived. 
One of the star from our sample, IRAS 12560, was previou sly ob- 
served in the CO J = 2— > 1 line bv lGroenewegen et alj JT997J). They 
derived a total mass-loss rate of 0.8x lO _6 M0yr _1 , which agrees 
fairly well with our estimate of 1.3x lO _6 M0yr _1 . 

Knowing the total mass-loss rate and the DPR allows us to 
estimate the gas-to-dust mass ratio ip. Assuming that the outflow 
velocities of the dust and gas are the same, ij) is simply the ratio of 
the total mass-loss rate to the DPR. Table|4] presents our estimates 
for tp using the DPRs from the radiative transfer models (ipo), and 
infrared colors (V>c). 



5 DISCUSSION 

5.1 Expansion velocities 

Lagad ecTet alj d2010h showed that the Halo carbon stars they ob- 
served have a rather high total mass-loss rate but a low expan- 
sion velocity. This result agrees with theoretical predictions by 
(Wachter et al. 2008) that the expansion velocity of carbon stars de- 
creases at lower metallicity while the mass-loss rate remains stable. 
They argue that these trends arise from less efficient dust forma- 
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Table 3. Dust pro duction rates (D PRs) and results from radiative transfer models for the observed stars and some comparison stars. For the Halo stars, 
these data are from | Lagadec et alj 1201(f) , except the last two entries, which are in the Halo betwee n SagDEG and the Galactic Center. Those data are from 
lLagadec et ailtoOar The modelling for Lynga 7 V 1 is new, based on data presented by Sloan et alj bOlfj) . 



Target 


Luminosity 


DPJ?DUSTY 


T e// 


Tin 


SiC/AmC 


t (1 1/xm) 


DPR col 


Distance 




(10 3 L o ) 


(10- 8 M o yr- 1 ) 


(K) 


(K) 




X10- 1 


(lO- 8 M yr- 1 ) 


kpc 


IRAS 04188 


8.0 


0.8 


2650 


1700 


0.10 


2.3 


1.0 


6.5 


IRAS 08427 


7.5 


2.3 


2800 


1200 


0.13 


2.5 


2.2 


5.5 


IRAS 1 1308 


16.8 


2.9 


2800 


1200 


0.00 


1.6 


2.4 


2.5 


IRAS 12560 


7.5 


0.3 


2800 


1200 


0.00 


1.1 


0.2 


12.0 


Lynga 7 VI 


9.0 


3.3 


2800 


1200 


0.15 


2.6 


2.1 


7.3 


IRAS 16339 


13.9 


1.2 


2800 


1400 


0.14 


2.3 


1.8 


4.9 


IRAS 18120 


8.0 


1.3 


2800 


1100 


0.15 


2.6 


1.3 


6.7 


IRAS 18384 


11.5 


2.6 


2800 


1200 


0.10 


1.1 


0.9 


14.5 


IRAS 19074 


13.3 


3.8 


2800 


1200 


0.20 


2.1 


1.3 


15.0 



Table 4. CO J=3-2 lines properties of the observed Halo stars, gas mass-loss rates and gas-to-dust mass ratio determined using radiative transfer models and 
infrared colors. 



Target v cxp Ico M gas 4>D i>c 

kms -1 Kkms" 1 lO" 6 M yr" 1 



IRAS 04188 


11.5 


0.77 


6.8 


846 


677 


IRAS 08427 


16.5 


0.55 


4.2 


184 


193 


IRAS 11308 


11.5 


10.86 


13 


456 


551 


IRAS 12560 


3.0 


0.14 


2.6 


866 


1300 


IRAS 16339 


8.5 


2.59 


11 


905 


603 


IRAS 18120 


6.5 


0.72 


5.4 


415 


415 



tion at low metallicity. As the dust accelerates the winds less effi- 
ciently at low metallicity, the expansion velocities of the envelopes 
are lower. But only a moderate amount of dust is needed to acceler- 
ate the wind. Once this threshold is reached, the mass-loss rate will 
remain constant. Fig. [4] shows the expansion velocity as a function 
of the total mass-loss rates for Galactic disc stars (open symbols) 
and our sample (filled symbols). T he values for stars in the Galac- 
tic disc are from lLoup etai] Jl993l) . The observations used for these 
work were also made in CO. 

At a given mass-loss rate, stars from the Halo and the Sagittar- 
ius Stream have a low expansion velocity. The star with an expan- 
sion velocity of 16.5 kms" 1 (IRAS 08427) belongs to the thick 
disc, as well as the two stars with expansion velocities of 11.5 
kms" 1 (IRAS 04188 and IRAS 11308). These stars do not appear 
as outliers in this figure, even though their expansion velocities are 
in the lower part of the observed range. 

The stars from the Sagittarius Stream and the Halo are clear 
outliers and have low expansion velocities for their mass-loss rates. 
This w ork thus extends the conclusion drawn by lLagadec et al.l 
(2009), who studied the DPRs only. The millimeter CO studies give 
a picture that somewhat contradicts the conclusions drawn from the 
infrared. While the infrared spectral properties are consistent with 
a metal-rich character, the low outflow velocities are much more 
typical of metal-poor stars. 



5.2 Gas-to-dust mass ratios 

Observations of CO with JCMT allow us to measure the total mass- 
loss rates of stars, while infrared observations enable us to deter- 



mine their DPRs. We can thus estimate their gas-to-dust mass ra- 
tio ip. As the metallicity of the individual stars we observed is not 
known we assume a solar value for the CO abundances, so that the 
total mass-loss rates we derived have to be seen as lower limits. 
Thus, the gas-to-dust mass ratios are also lower limits. The esti- 
mated values of ip range from ~200 (the value typically assumed 
for Galactic stars) to 1300. The upper range of our estimates for ip 
is high, as might be expected for metal-poor stars, where we ex- 
pect less dust to form for a given amount of gas. IRAS 08427 has 
a gas-to-dust ratio close to the canonical Galactic value of ~200. 
This star belongs to the thick disc and has the highest expansion 
velocity and the lowest ip in our sample. 

Some caution is warranted when drawing conclusions from 
such small samples, but these observations are pointing to high gas- 
to-dust mass ratios and low expansion velocities in carbon stars in 
the Halo. To generalize, the gas-to-dust mass ratios of carbon stars 
appear to increase at lower metalli cities. Such a trend is expected 
from a theoretical point of view dvanLoonl hoOO; Wacht er et alj 
2008), as dust formation is expected to be less efficient in metal- 
poor environments. 

IRAS 04188 is an outlier, as it is a thick disc star with a high 
gas-to-dust ratio. It is interesting to note that the spectral energy dis- 
tribution of this star could not be fitted with the same standard pa- 
rameters as the other stars. A rather cool star (2650 K) was needed 
for the fit, as well as very hot dust (1700 K) for the dust inner ra- 
dius. 

For some of the thick disc stars we observed, ip is quite large. 
We measured ip by assuming that the drift velocity (the difference 
between the dust and the gas velocity) is equal to zero. If the out- 
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Figure 4. Expansion velocity (in krns" 1 ) as a function of the total mass- 
loss rate (in MQyr" 1 ). Open symbols are carbon stars in the Galactic disc. 
Filled symbols are the stars from our sample, with squares for the thick disc 
and circles for the Halo. 

flow velocities of the gas and dust were decoupled, then the actual 
value for ip would be smaller than indicated here. Such a decou- 
pling would require less efficient transfer of momentum from the 
radiatively accelerated dust grains to the gas in more metal-poor 
environments, which might be possible if the grain-size distribution 
were shifted to larger grains. We note that the coupling efficiency 
is determined primarily by the gas density. 

5.3 Metallicity of the Galactic halo carbon stars 

Fig.|5]displays the equivalent width of the 7.5 /im absorption band 
from C 2H2 as a function of the strength of the SiC dust emission 
feature. lLagadec et al.l J2009h introduced this plot and showed that 
it can diagnose the metallicity of carbon stars observed spectro- 
scopically in the infrared. Despite some scatter and mixing between 
the populations, the carbon stars from the Milky Way are clearly 
concentrated toward the lower right corner. Moving progressively 
up and to the left, one moves through the region dominated by the 
LMC, followed by the SMC. The Fornax stars trace the same region 
as the SMC, consistent with the similar metallicities of the carbon 
stars sample d with the IRS to t hose in the Magellanic Clouds, as 
described bv lSloan et alj d2012h . The position of stars on this dia- 
gram depends also on its initial mass, sp that it can be use only as 
an indication of the metallicity in the present case. 

The carbon stars from our sample, both in the thick disc and 
Halo, as well as our two comparison stars in the foreground of the 
Sgr dSph and the one star near Lynga 7 VI, occupy the same region 
as the Galactic carbon stars. 

As mentioned in Section[3] decreasing the initial metallicity 
tends to lead to a higher C/O ratio and thus stronger C2H2. This 
would also decrease the initial abundance of Si and thus decrease 
the SiC strength. But silicon, being a heavy a-element, it is not 
synthetized in AGB stars. The location of the star on the diagram 
Fig.|5]thus not only depends on [Fe/H] but also on the initial a- 
element abundance. 



Figure 5. The equivalent width of the C2H2 absorption band at 7.5 fim 
as a function of the ratio of the strength of the SiC dust emission feature 
at 11.3 /im compared to the underlying continuum. The symbols are as 
defined before, with the largest open circles representing the Halo stars, 
and the intermediate open circles representing stars in the thick disc. 

On the other hand, observations of C-rich post- AGB stars with 
the 21 microns feature in the Galaxy and the Magellanic Clouds 
(Decin et al., 1998; De Smedt et al., 2012 and references therein) in- 
dicate that metal-poor carbon stars can have significantly enhanced 
C, O and Si abundances. Such an enhancement would explain the 
strong SiC feature observed in the metal poor carbon stars we ob- 
served in the halo. However, it is not clear why metal-poor carbon 
stars would produce more Si in the halo than in other galaxies of 
the Local Group, as we have shown AGB stars in the SMC, LMC, 
Fornax and Sculptor had a weaker SiC feature than Galactic stars. 

Enrichment of the ISM by Type II supernovae, which are 
known producers of a elements is a tempting explanation, but the 
known behav i or of [ q/Fe] as a function of [Fe/H], as illustrated by 
iTolstov et alj j2009L see their Fig. 11) is not consistent with our 
spectra. The [a/Fe] abundance increases by ~one decade while 
[Fe/H] drops by two decades. In terms of absolute abundance, 
[a/H] is decreasing at lower metallicities. The bottom panels of 
Fig. 2 contradict this expected behavior. They show that the SiC 
features in our Halo and thick disc spectra are strong compared to 
most of the other spectra in absolute terms and not just as a ratio to 
amorphous carbon, as shown in the upper panels. 

The two Halo stars observed with Spitzer and ICMT appear 
to have a low expansion velocity (an indication that they are metal- 
poor) and a very strong SiC feature. Those stars could thus either 
be similar to Galactic disc stars or metal-poor and a-element en- 
hanced. 

5.4 Origin of the Galactic halo carbon stars 

The typical age of halo stars is 10 Gyr, so no AGB stars more mas- 
sive than 1 Mq formed in the Halo should be observable now. Mod- 
els of AGB stars with such a low mass do not dredge up enough 
carbon to make the stars carbon-rich. But carbon-rich stars are ob- 
served in the halo (Totten & Irwin, 1998; Mauron, 2008). The fact 
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that we observe carbon stars in the halo could be a clue that these 
models are not right, or that these stars were formed outside of 
the halo. Comparing the luminosity of the star we observed with 
Fig. 20 from (Vassiliadis & Woods 1993) indicates that the initial 
masses of the stars in our sample are between 2 and 3M0 . This is 
one more indication that those stars might have formed outside the 
Galactic halo. We propose two possible external origins for these 
halo carbon stars. 

It has been proposed that carbon AGB stars in the halo could 
originate from the tidal disruption of Sgr dSph Galaxy, orbiting the 
Milky Way with a period of about 1 Gyr (Ibata et al. 2001). Their 
model of the Sgr dSph stream shows a good agreement between the 
orbit of the stream and those of observed carbon stars. 

Out of the three halo stars we observed with the JCMT, one, 
IRAS 12560, likely belongs to the Sagittarius stream (Lagadec et 
al. 2010) and thus originates from Sgr dSph. These three Halo 
stars have a very low envelope expansion velocity, indicative of 
a metal-poor nature (Lagadec et al. 2010). The Spitzer spectra we 
present here show that those stars present strong C2H2 absorption 
bands, as observed in metal-poor carbon stars (see e.g. Sloan et 
al. 2009). These two arguments strongly support the hypothesis 
that these stars are metal-poor. As shown in Section [531 the halo 
carbon stars we observed could be enhanced in a-elements. They 
would thus have formed in a burst-like stellar population. As they 
are carbon stars, this population is necessarly younger than the halo 
itself. The Sgr dSph have such a stellar population. The halo car- 
bon stars we observed may thus come from the Sgr dSph galaxy 
or a similar galaxy orbiting the Milky Way. The fact that the dust 
properties of the Sgr dSph carbon stars are very similar (Lagadec et 
al. 2009) supports this argument. This galaxy contains a metal-rich 
population (Kniazev et al. 2007) but it is not known whether this 
population is represented in the tidal tails. 

But those stars could also be similar to Galactic disc stars 
(Section |5"31 >, This leads naturally to the question of how stars 
formed near the Galactic plane could wind up several kpc above 
or below it. They could be runaways from the thin disc. One inter- 
esting scenario is that they may have been ejected from the break- 
up of a binary, possibly as a result of the interaction of a binary 
and a single star, as suggested for more massive run aways (e.g., 
iHoogerwerf et al.ll200ll ; lFuiii & Portegies Zwartll201 ll) . While car- 
bon stars are bright and easily detected, they represent perhaps one 
thousandth of a star's total lifetime. This would point to a larger 
population of less evolved intermediate-mass stars now in the Halo, 
but relatively metal-rich. 

To discriminate the two possibilities, we would need to obtain 
measurements of [Fe/H] and different a-elements abundances for 
the stars we observed. Such measurements are made difficult by the 
fact that those stars are very faint in the optical. 

A final possibilty is that these are extrinsic carbon stars, with 
the carbon enrichment caused by mass transfer from a more mas- 
sive, faster evolving companion. 



6 CONCLUSIONS 

We presented Spitzer IRS spectra of four Galactic disc and halo 
carbon stars. Typical features of carbon stars with SiC dust emis- 
sion and C2H2 molecular absorption are seen in the four spectra. 
Dust radiative transfer models enabled us to determine the dust- 
production rates for these stars. We measured the gas mass-loss 
rates and expansion velocities from the stars via modeling of the 
CO J = 3 — y 2 transition. The gas-to-dust mass ratio ip was then 



measured by simply dividing the gas and dust mass-loss rates. The 
measured mass-loss rates and tp are in the range 1-4 lO~ 6 M0yr -1 
and 200-1300 respectively. The stars in the halo have a low expan- 
sion velocity (around 7 km/s), strong C2H2 molecular absorption 
bands, strong SiC emission and a high gas-to-dust mass ratios. 

The strong SiC emission we observe could be due to an over- 
abundance of a-elements (like silicon) in the Halo, that these stars 
produce silicon or to the fact that those stars metallicities similar 
to the Galactic disc. In the latter case, they may have escaped from 
the disc. The low expansion velocities and strong C2H2 absorption 
bands are arguments in favor of the former. However, explaining 
the a enhancement in the Halo is a challenge. We conclude that 
these spectra present us with some intriguing contradictions, with 
the dust properties being similar to metal-rich carbon stars and the 
gas properties similar to metal-poor carbon stars. 

Those AGB stars could have formed outside of the Halo, as 
the mass of these carbon stars is inconsistent with a formation in 
the Halo. The halo carbon AGB stars we observed could thus have 
two possible origins. They could have formed in a galaxy similar to 
Sgr dSph, orbiting the Milky Way and being gradually stripped by 
tidal forces or they could have escaped from the Galactic disc. 

This work has shown the necessity to combine mid-infrared 
spectroscopy with millimeter CO observations to fully characterize 
the dust production and gas mass-loss from AGB stars. Observa- 
tions of metal-poor AGB stars in Local Group dwarf galaxies are 
necessary to study the dependence of this mass-loss on metallicity. 
The Atacama Large Millimeter/submillimeter Array (ALMA) will 
allow such observations. Finally near-infrared spectroscopy of the 
individual stars will be necessary to precisely determine the metal- 
licity and a-element abundances of the observed stars. 
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